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In recent years, many variants of percolation have been used to study network structure and the 
behavior of processes spreading on networks. These include bond percolation, site percolation, k- 
core percolation, bootstrap percolation, the generalized epidemic process, and the Watts Threshold 
Model (WTM). We show that — except for bond percolation — each of these processes arises 
as a special case of the WTM and bond percolation arises from a small modification. In fact 
“heterogeneous k- core percolation”, a corresponding “heterogeneous bootstrap percolation” model, 
and the generalized epidemic process are completely equivalent to one another and the WTM. We 
further show that a natural generalization of the WTM in which individuals “transmit” or “send a 
message” to their neighbors with some probability less than 1 can be reformulated in terms of the 
WTM, and so this apparent generalization is in fact not more general. Finally, we show that in 
bond percolation, finding the set of nodes in the component containing a given node is equivalent 
to finding the set of nodes activated if that node is initially activated and the node thresholds 
are chosen from the appropriate distribution. A consequence of these results is that mathematical 
techniques developed for the WTM apply to these other models as well, and techniques that were 
developed for some particular case may in fact apply much more generally. 


INTRODUCTION 

To understand processes spreading on a static network 
G, researchers frequently investigate how G behaves un¬ 
der percolation. Percolation comes in many flavors, and 
the information we gain depends on which variety we 
choose. Most frequently, we study bond or site perco¬ 
lation, but researchers have also found that k -core per¬ 
colation, bootstrap percolation, the generalized epidemic 
process, and the Watts threshold model (WTM) provide 
valuable insights m m uni usi [in m]. These processes 
are closely related, and indeed similar mathematical ap¬ 
proaches have been used to study several of these pro¬ 
cesses mm- Our main result is that all of these (and 
some related) processes can be derived as special cases 
of the WTM, and in fact several of these are completely 
equivalent to the WTM. 

Much of the motivation for studying percolation pro¬ 
cesses comes from trying to understand spreading pro¬ 
cesses in networks. If we consider systems in which nodes 
change status in response to the status of their neigh¬ 
bors, and the potential path of statuses they can have 
is acyclic (that is they can never return to a previous 
status), then many variants of percolation can be ap¬ 
plied. This is commonly used for SIR disease, in which 
an individual can be infected by an infected neighbor. 
However, much recent work has focused on the spread 
of “social contagion” or “complex contagions” [121 i T3] in 
which multiple transmissions may be required in order to 
cause “infection”. Sometimes this is presented as assign¬ 
ing each node a threshold r u such that u becomes infected 
once r u neighbors are “infected”. Other times this is 
presented as a reduction (or increase) in the probability 
that a neighbor will transmit as an individual encoun¬ 
ters more infected individuals. This models the idea that 
after hearing seemingly independent “confirmation” of a 


rumor people may be more likely to believe and spread it, 
or after seeing multiple people engaging buying a prod¬ 
uct someone is more likely to perceive a consensus and 
buy the product as well. Some experimental evidence of 
this has been found 1 lj j33j. 

We briefly review the processes we will study: In bond 
percolation, some edges are independently selected with 
uniform probability p to be retained while the remaining 
edges are deleted (with probability 1 — p). Similarly in 
site percolation, some nodes are randomly selected with 
probability p and the remaining nodes are deleted. Typ¬ 
ically our interest is in identifying the nodes in the con¬ 
nected components of the residual network, and whether 
a “giant” component exists (that is a component whose 
size is proportional to the network size in the infinite 
network limit). 

Bond percolation and site percolation often show up 
in the study of SIR disease spread where a single trans¬ 
mission suffices to cause infection isi nni I 2 H 123 na eu- 
[H El 131131133- There is an exact equivalence between 
the spread of an SIR disese and bond percolation, and so 
much has been learned about the threshold, scaling prop¬ 
erties, and dynamics of an SIR disease by studying the 
corresponding percolation odel. This percolation equiv¬ 
alence is based on the fact that an edge either exists or 
does not in percolation, while in disease spread if the 
edge transmits, the receiving node becomes infected. 

In k-core percolation all nodes with degree less than 
some specified k are removed. This removal may reduce 
some nodes’ degrees below k. If so, these are removed. 
This “pruning” process repeats until reaching a state in 
which all nodes have degree at least k. This remaining 
network is called the “fc-core” of the network. It is seen 
to have hybrid phase transitions, with a square root-type 
scaling on one side of a transition followed by a discontin¬ 
uous jump HE]- In a variant, “heterogeneous fc-core” 


2 


percolation [7J, each node is assigned its own threshold 
value and deleted if its degree goes below the threshold. 
We note that many authors have used the term “boot¬ 
strap percolation” to denote k- core percolation, and in¬ 
deed this appears to be the original term [2JQJ3HD], but 
we reserve “bootstrap percolation” for a closely related 
dual process. The k- core has been applied to many prob¬ 
lems, including understanding failure of a physical system 
under strain |38| , network visualization [5] , identification 
of the component of a network responsible for establish¬ 
ing a disease |24| . and more generally for understanding 
the structure of a network m- 

In bootstrap percolation (introduced in |2j where it is 
called “diffusion percolation”) a collection of nodes is ini¬ 
tially “activated”. Then any inactive node with at least 
to active neighbors becomes active. The process repeats 
until all remaining inactive nodes have fewer than to ac¬ 
tive neighbors. It was initially introduced to model the 
spread of a water-filled crack in a rock. It has received 
considerable study on lattices miss, and its behavior 
in large random networks has been the subject of some 
more recent analysis [B]. Like k- core percolation it is 
seen to have a hybrid phase transition. We introduce 
a natural generalization analogous to heterogeneous k- 
core percolation in which each node is assigned its own 
threshold. This “heterogeneous bootstrap percolation” 
does not appear to have been studied previously. 

In the generalized epidemic process (GEP) [5] 1161 122| . 
we think of an infection spreading through the network. 
If a node has a single infected neighbor, its probability 
of becoming infected is pi. If it escapes infection, but 
a second neighbor becomes infected, then its probability 
of becoming infected is P 2 . This repeats and the proba¬ 
bility of successful transmission on the m-th neighbor’s 
infection is p m . If p m = p for all to, then this is the net¬ 
work version of the classical Reed-Frost model [T] for a 
susceptible-infected-recovered disease [5] ■ If Pm decreases 
as to increases, this could model decreasing susceptibility 
due to an improved immune response as exposures accu¬ 
mulate, or it could simply represent pre-existing hetero¬ 
geneities in susceptibility that are revealed as the number 
of exposures increases. An increasing p m would model 
some synergistic or cumulative effect of exposures as seen 
in “complex contagions” jT2]. For comparison with other 
models, we allow p m to depend on d u , the degree of node 
u. 

In the WTM 17, [39], each node u is assigned an in¬ 
dividual threshold r u which we assume is assigned to u 
independently at random, with a probability that may 
depend on its degree d u . The probability node u has 
a given r is given by P(r u = r\d u ) = q{r\d u ). A node 
begins either active or inactive. If an inactive node u 
has at least r u active neighbors then it becomes active. 
We assume that the initially active nodes may be cho¬ 
sen independently at random (which can be modelled by 
having q(r\d) > 0 for some r < 0), or they may be cho¬ 


sen by some other rule, in which case we treat the set of 
initially active nodes as an input to the algorithm. Often 
a common threshold r* is chosen so P(r u = r*) = 1 or a 
common fraction p* is chosen so P{r u = [p*d u ] |d u ) = 1. 
As described above, this is frequently used to model so¬ 
cial contagions. In [35] it was conjectured that for a 
global cascade to occur from an infinitesimally small ini¬ 
tial proportion active, a giant component of nodes with 
r = 1 would need to exist. This is true in random Con¬ 
figuration Model networks, but false in random clustered 
networks [5111 . As with bootstrap and k -core percolation, 
this is known to exhibit hybrid bifurcations [13- 

In these generalized percolation processes, typically we 
are interested in the final set of active nodes, but some¬ 
times we may be interested in the temporal dynamics as 
these nodes become active 0(30]. If we are interested in 
the temporal dynamics, then we must assign additional 
rules for how long it takes for a node to become active. 
Although the timing will depend on the details of the 
additional rules, the final set of active nodes is uniquely 
determined once the network, thresholds, and initially 
active nodes are chosen. For our purposes we focus just 
on the final state. 

We will show that by appropriately choosing the dis¬ 
tribution of r and the initial set of active nodes, we can 
recover other versions of percolation from the WTM, in¬ 
cluding site percolation, k -core percolation, bootstrap 
percolation and the GEP. Going a step further, we show 
that the heterogeneous fc-core of a network, the deleted 
nodes in heterogeneous bootstrap percolation, and the 
set of “infected nodes” in the GEP are in fact all equiva¬ 
lent to the set of active nodes emerging from the WTM. 
That is, given one model and the corresponding distri¬ 
bution of thresholds, we can define the distribution of 
thresholds of the other models to yield the same sets of 
nodes with the same probabilities. A natural generaliza¬ 
tion of the WTM to consider has each node “transmit¬ 
ting” or “passing a message” with some fixed probability 
T. We show that by modifying the threshold distribution 
the original WTM (with T = 1) can recover the same 
outcomes as for any other To < 1, and thus allowing for 
T < 1 does not enlarge the set of possible outcomes. 

Finally, we investigate the relation with bond percola¬ 
tion. If our interest in bond percolation is to identify the 
connected component containing a given node rt, then we 
can find this component using the WTM with u as the 
initially active node and appropriate threshold distribu¬ 
tion. To find all connected components, we can start the 
WTM with one initially active node, run it to completion, 
and then choose a remaining inactive node and rerun the 
WTM, iterating until no inactive nodes remain. The set 
of nodes that are activated in each pass correspond ex¬ 
actly to the components found in bond percolation. 
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ANALYSIS 

We begin by explicitly describing an algorithm which 
implements the WTM. Each node is assigned a weight w 
uniformly between 0 and 1. We use a model-dependent 
function dist func to convert the weight into a threshold 
r, possibly depending on the degree of the node. Typi¬ 
cally we choose the function to return the largest value 
r u such that JJr=-oc ?( r M) < If there are specified 
initially active nodes, they are given a threshold of zero. 
Alternately we can allow the randomly assigned thresh¬ 
old to permit values r u < 0 in which case these nodes 
are initially active, and the iterative process begins. For 
each active node, we reduce the threshold of any inactive 
neighbor by 1. If a node’s threshold reaches 0 it activates. 
Pseudocode for the algorithm is given in the appendix in 
figure [6] 

Once the random thresholds and index nodes are set, 
the final outcome of the WTM is deterministic. To show 
that the other percolation processes give the same behav¬ 
ior, we will show how to structure these processes to start 
from the same random weights w u and deterministically 
yield a final state that is identical to the state found by 
the WTM for some threshold distribution. 

Site Percolation 

In site percolation, each node is retained with prob¬ 
ability p or deleted with probability 1 — p. To simu¬ 
late site percolation, we can generate a random number 
w u £ (0,1) independently and uniformly at random for 
each node u. If w u < p (which occurs with probability p) 
we keep u , otherwise we delete it. It is straightforward 
to see that this is identical to the algorithm in Figure [6] 
if the threshold is set to be r u = 0 whenever w u < p and 
r u = d + 1 otherwise. In this case, with probability p the 
node has threshold 0 and so is initially active, while with 
probability 1 — p it has threshold d + 1, and so can never 
become active as it will have at most d active neighbors. 
Thus, nodes are retained in site percolation iff they are 
active in the WTM. This is demonstrated in figure [T] 

k -core Percolation 

We now consider k- core and heterogeneous k- core per¬ 
colation. The classical fc-core percolation is determinis¬ 
tic: each node with fewer than k neighbors is deleted. 
This iterates until all remaining nodes have at least k 
neighbors among the remaining nodes. To reproduce this 
with the WTM, we set r u = d u — k + 1 regardless of w u . 

With this threshold, all nodes with d u < k activate im¬ 
mediately in the WTM. In k- core percolation these same 
nodes are immediately deleted. For a given node u not 
in this set, let the number of neighbors activated/deleted 



FIG. 1: Comparison of site percolation and WTM on honey¬ 
comb lattice. Results for p = 0.4 on left and p = 0.8 on right, 
(top) Each node is assigned a weight, (middle) Site percola¬ 
tion: If the weight is less than p , the node is kept, otherwise 
it is deleted, (bottom) WTM: if the weight is less than p it 
is given a threshold of 0. Otherwise it is given d + 1. Those 
with threshold 0 are shown in color, and activate immediately. 
Those with threshold larger than their degree are uncolored 
and never activate. 

be denoted n u . In the WTM, any remaining node with 
d u — k < n u then activates. In k- core percolation, any 
node with d u — n u < k is deleted. Again, these nodes 
are the same. Iterating as shown in figure [2j the set of 
activated nodes in the WTM is the set of deleted nodes 
in k- core percolation. 

We can repeat this for heterogeneous k-core percola¬ 
tion. We assign weights w u to each node and map that 
to a heterogeneous k- core threshold k u . We can map this 
weight to a WTM threshold such that if the node is as¬ 
signed a given k u , it is assigned r u = d u — k u + 1 for the 
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FIG. 2: Comparison of the first steps of k- core percolation and 
the WTM for the karate club graph [30]. (Left) k- core perco¬ 
lation with k = 4. (top) The original network, (middle) The 
first step of k- core percolation, (bottom) The second step. 
(Right) The WTM (top) Thresholds of d u — 4 + 1. (middle) 
Nodes with thresholds r u < 0 are activated, (bottom) The 
second step. At each step, the activated nodes of the WTM 
are exactly the deleted nodes in k- core percolation. 

WTM. Then the WTM and heterogeneous k- core perco¬ 
lation are equivalent: a node is deleted in heterogeneous 
k- core percolation iff it is activated in the WTM. 



FIG. 3: A comparison of heterogeneous bootstrap and hetero¬ 
geneous k- core percolation for the social network of dolphins 
observed by [26]. (Left) Heterogeneous bootstrap percolation: 
(top) thresholds for activation, [d„/3j. (middle) first step: all 
nodes of degree 1 or 2 are activated, (bottom) second step: 
nodes that now reach their threshold are activated. (Right) 
Heterogeneous k- core percolation: (top) thresholds for dele¬ 
tion, d u — [d u / 3J + 1 . (middle) first step: all nodes of degree 
1 or 2 are deleted, (bottom) second step: nodes that now 
reach their threshold are deleted. The nodes deleted at each 
stage of k- core percolation correspond exactly to the nodes 
activated at the same stage of bootstrap percolation. 


Bootstrap Percolation 

In bootstrap percolation, some initial nodes are acti¬ 
vated, and nodes become active once they have at least 
k active neighbors (k is the same for all nodes). This 
is similar to k- core percolation, but k- core percolation is 
subtractive while bootstrap percolation is additive 00. 

We consider bootstrap percolation with a set Iq of ini¬ 
tially active nodes, and compare it to the WTM with 
r u = k for all nodes except the nodes in Iq which are 
initially active. Following a similar argument to the 
WTM/fc-core percolation equivalence, we see that with 
this definition, the WTM adds nodes to the system ex¬ 
actly when bootstrap percolation does. 

If we consider heterogeneous bootstrap percolation, 
then a similar argument also shows that it is equivalent 
to the WTM. Because of the correspondence between the 
WTM and heterogeneous fc-core percolation, this means 
that heterogeneous bootstrap percolation is equivalent to 


heterogeneous k- core percolation, with the deleted nodes 
in heterogeneous k- core percolation matching the acti¬ 
vated nodes in bootstrap percolation. 

At first glance, this contrasts with observations of 0. 
They showed that the fc-core and the activated nodes in 
bootstrap percolation are not the same and can have dif¬ 
ferent internal structure. In fact, the distinction between 
the two turns out to be that the nodes defined to be 
active for the bootstrap version are the nodes deleted 
in the k- core version. They are complementary pro¬ 
cesses. Any behavior observed in heterogeneous k- core 
percolation can be observed in heterogeneous bootstrap 
percolation the inactivated nodes of heterogeneous boot¬ 
strap percolation, while any behavior observed in the ac¬ 
tivated nodes of bootstrap percolation can be found in 
the deleted nodes of fc-core percolation. This equivalence 
is previously known J2J. 

Figure [3] demonstrates the equivalence between hetero¬ 
geneous bootstrap and heterogeneous k- core percolation. 
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Generalized Epidemic Process 

We now consider the generalized epidemic process 
(GEP) [9j 122, for which the m-th “infected” neighbor 
infects node u (given that the previous m — 1 did not) 
with probability p m (d u ). Our approach resembles the 
“Sellke construction” [3B] of a simple epidemic model in 
a fully-mixed population. In a standard fully-mixed epi¬ 
demic simulation, an individual that is susceptible at the 
start of a short time interval becomes infected with a 
probability proportional to the number of infected indi¬ 
viduals. In the Sellke construction formulation, however, 
we assume we know in advance for each individual the 
cumulative amount of exposure it will receive before be¬ 
coming infected (this is a random number chosen from 
an exponential distribution). We then begin the spread 
with some initial infections and when (or if) the exposure 
reaches that threshold the individual becomes infected. 

We will now study the network-based GEP by a sim¬ 
ilar approach. The probability the first m — 1 in¬ 
fected neighbors do not infect u but the m-th does is 
Pm{du)Yl2=l^ — Pm(d u ))- We simply assign a random 
number w u G (0,1) and map this to m u . Thus for any 
given node, it will become infected upon the infection 
of its m-th neighbor with probability p m {d u ) IIa=i(1 ~ 
Pm(d u )) independently of other nodes and independently 
of whether we will calculate m u in advance or simply 
accept or reject infection with probability Pm{d u ) as it 
accumulates infected neighbors. 

For the WTM we use the same mapping from w u to r u 
so r u = m u . The node u activates exactly after the ru¬ 
th neighbor activates, while in the GEP u is infected at 
exactly the same step. Thus any GEP can be expressed 
as a WTM. Showing the inverse is straightforward, and 
so the GEP and WTM are equivalent. If we do not allow 
p m to depend on d u (as in the original version), then this 
is a special case of the WTM. 

Bond Percolation 

We finally consider Bond Percolation. Typically in 
bond percolation, we can consider the edges in any order, 
choosing to keep each edge with probability p or delete 
it with probability 1—p independently of the others. We 
then identify the connected components of the network. 

We will focus our attention just on identifying which 
nodes form connected components after bond percola¬ 
tion; we are not interested in which edges exist within the 
components. In figure [4] we compare a bond percolation 
approach to finding the component containing a partic¬ 
ular node with a WTM approach for finding the same 
component. We first perform bond percolation. We then 
select an initial node (highlighted in the figure), and fol¬ 
low edges out from that node in the percolated network 
to find its component. Nodes are labeled with r where r 



FIG. 4: (top left) Original network, with initial node high¬ 
lighted. (top right) Percolated network with component of 
initial node highlighted taking p = 0.51. Edges shown within 
a component are red, while edges in other components are 
black, (bottom left) WTM outcome with thresholds found 
from the percolated network using a breadth-first-search for 
the WTM. (bottom right) WTM outcome with thresholds 
from the percolated network using a depth-first search for the 
WTM. In both WTM plots, the edges that were responsible 
for the activation of a node are shown in red. Edges which 
were never considered are shown dashed in black. 

is the number of edges of the original network that were 
encountered (but deleted) prior to an undeleted edge. 

We can think of this as being indistinguishable from 
selecting an initial node, following edges out from that 
node in some order, where each time an edge is consid¬ 
ered, it is deleted with probability 1 — p or followed with 
probability p. The probability that the first r edges to a 
node are deleted but the next is not is p[\ — p) r . 

We compare this with the WTM with a threshold of 
t u = r u + 1. The activated nodes are identical to the 
component found using bond percolation. In general, 
assigning nodes a threshold of r where r > 1 is taken 
with probability p{ 1 — p) T will yield a set of active 
nodes from an initially active node which come from the 
same distribution as the component of that node follow¬ 
ing bond percolation. 

In fact, we can generalize this approach to find all the 
components. The steps in our process are to begin with 
a network, and assign thresholds using a geometric dis¬ 
tribution: for a threshold of r the probability of r is 
p( 1 — p) T -1 . We then select a node and successively add 
nodes to its component once their threshold number of 
neighbors have been visited. This process is likely to 
terminate without exploring all nodes. If this happens, 
we iteratively select a new node and add nodes to its 
component whenever their threshold number nodes have 
been visited (either in this stage or while building a pre¬ 
vious component). The resulting components match the 
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components observed in bond percolation. Nodes are ac¬ 
tivated exactly when they are added to a component in 
the bond percolation, and identifying in which iteration 
they are activated tells us which component they are part 
of. The breadth-first search figure is the implementation 
of the WTM shown in Figure [6j but we highlight that an 
alternate implementation with a depth-first search would 
yield the same outcomes. As long as the initial nodes 
of each pass are chosen in the same order, the precise 
details of the search algorithm do not determine which 
additional nodes belong to the component identified in 
that pass. 




FIG. 5: Activated clusters found using depth-first (left) and 
breadth-first (right) searching using the WTM with a thresh¬ 
old of r occuring with probability p( 1 —p) T ~ 1 (independently 
of d) and p = 0.51 for a 9 x 9 lattice. The number at each 
node is its threshold. The circled nodes are the initial nodes 
chosen for each cluster. The bottom left node is chosen first, 
and its cluster traced out. The next cluster is initialized by 
the bottom-most of the left-most remaining nodes. Thick col¬ 
ored edges formed the final interaction that caused activation. 
Non-existent edges failed to cause activation (but moved the 
node closer to its threshold). Dashed black edges were not 
tested because both nodes were already active when the edge 
was considered. The clusters remain the same for both search 
orders (but edges change). 

To arrive at bond percolation, the thresholds for the 
WTM process are assigned from a geometric distribu¬ 
tion. It would be interesting to study whether a different 
distribution could be interpreted in the context of a gen¬ 
eralized bond percolation. 

DISCUSSION 

Many percolation processes have been studied in net¬ 
works. We have shown that site percolation, bootstrap 
percolation, k- core percolation, and the GEP are all spe¬ 
cial cases of the WTM. In fact, the GEP we consider is 
equivalent to the WTM, and if we allow a node-specific 
threshold then both bootstrap and k- core percolation are 
also equivalent. Which one should be considered the 
“base” model is a matter of personal choice. 

Bond percolation is closely related to the WTM, but to 
arrive at an equivalent model, the WTM assigns thresh¬ 
olds from a geometric distribution, activates a node, fol¬ 
lows the WTM process to completion, and then activates 


another node. The successive sets of activated nodes oc¬ 
cur with the same probability as would be found in bond 
percolation. 

We have further shown that generalizing the WTM 
to allow for a homogeneous transmission probability T 
from active nodes to neighboring inactive nodes results 
in a model which can be thought of as a special case 
of the WTM. Thus the potential space of models is not 
increased by this modification. 

This commonality helps to explain why similar behav¬ 
iors are observed and similar mathematical methods ap¬ 
ply to these different processes. 
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Appendix: Algorithm 

In figure [6] we give pseudocode for the WTM algo¬ 
rithm. Other implementations are possible (this one is 
based on a breadth-first search, but for example, a depth- 
first search could also be used). The choice of the func¬ 
tion dist func which maps a randomly chosen weight from 
(0,1) to a threshold allows us to match other percolation 
models. 
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Input: Input network G, function generating numbers from 
a distribution dist func, and set of initially active nodes Jo. 
Output: Set ActivatedNodes of activated nodes. 


function WTMAssignVVeights(G) 

for u in G. nodes do 

Assign weightju] uniformly from (0,1) 
return weight 


function WTM_Assign_Thresh(G, distfunc, weight, Jo) 
for u in G.nodes do 
if u in Iq then 
thresh [it] <— 0 
else 

thresh[it] t— dist func(G.degree(u), weight[u]) 
return thresh 


function WTM_Process(G, thresh) 

CurrentNodes t— set of nodes in G with thresh < 0 
ActivatedNodes t— set of nodes in CurrentNodes 
while CurrentNodes is not empty do 
NextNodes <— emptySet 
for u in CurrentNodes do 
for v in G.neighbors(ti) do 
if v qL ActiveNodes then 
ActiveN odes. add (w) 
NextNodes.add(u) 

CurrentNodes t— NextNodes 
return ActiveNodes 


function WTM(G, dist func, 7o) 

weight = WTMAssignWeights(G) 
thresh = WTM Assign Thresh(G,dist func, weight, 7o) 
ActivatedNodes = WTM_Process(G, thresh) 
return ActivatedNodes 


FIG. 6: The steps of the WTM algorithm. The main al¬ 
gorithm is the final function given. The other functions are 
called by the main algorithm. First the weights are assigned 
randomly, and then the weights are mapped (deterministi¬ 
cally) to a threshold and algorithm proceeds iteratively (and 
deterministically). The appropriate choice of dist func allows 
us to select between the different models. 
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